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Abstract 

We consider a competing risks model, in which system failures are due to one out of 
two mutually exclusive causes, formulated within the framework of shock models driven 
by bivariate Poisson process. We obtain the failure densities and the survival functions as 
well as other related quantities under three different schemes. Namely, system failures are 
assumed to occur at the first instant in which a random constant threshold is reached by 
(a) the sum of received shocks, (b) the minimum of shocks, (c) the maximum of shocks. 

1 Introduction 

The classical competing risks model deals with failure times subject to multiple causes of 
failure. This model is of specific interest in various applied fields such as survival analysis and 
reliability theory. Indeed, it is appropriate for instance for describing failures of organisms or 
devices in the presence of different types of risks. Usually this model involves an observable 
pair (T, 5), where T is the time of failure and S describes cause or type of failure. General 
properties of this model can be found in the literature (see Bedford and Cooke [3J and Crowder 
[5]). For a list of references and recent results on competing risks model see Di Crescenzo and 
Longobardi [S], and other contributions in [TJ. 

In this paper we introduce a formulation of the competing risks model within the frame- 
work of stochastic shock models. These have been introduced and studied with the aim of 
describing systems subject to shocks occurring randomly such as events in counting processes. 
About 40 years ago A.W. Marshall and I. Olkin began to investigate probabilistic shock mod- 
els. In particular, in 25J and 26 they obtained certain bivariate exponential distributions 
from fatal and non-fatal shock models. Later, by means of the total positivity theory, Esary et 
al. [H] introduced a classical preservation problem in shock models, i.e. to show that certain 
properties of the discrete distribution of the number of shocks that leads a system to failure 
are reflected in corresponding properties of the continuous distribution of the system life- 
time. This investigation line has been followed by numerous authors including, for instance, 
El-Ncweihi et al. [TP] . We also recall the contributions by Gottlieb [TB] and by Ghosh and 
Ebrahimi [15] , concerning the increasing failure rate property of the life distribution in shock 
models. Furthermore, preservations results for partial stochastic orderings and for classes 
of failure distributions in various types of shock models have been obtained by Klefsjo [21"] . 
Fagiuoli and Pellerey [12], [13], [14], Kebir [19], Kochar [22], Pellerey [29], Singh and Jain [34]. 
Moreover, cumulative shock models in univariate and multivariate cases have been considered 
by Gut [IS], Kijima and Nakagawa (2U], Perez-Ocon and Gamiz-Perez [31] , [52"] . Shaked and 
Shanthikumar [33]. 

Recently, new general classes of shock models of interest in reliability theory and survival 
analysis have been proposed, in which failure is due to the competing causes of degradation and 
trauma (see Lehmann [23] , where the survival function of the failure time has been computed). 
This research direction is motivated by the circumstance that failure mechanisms can often 
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be ascribed to an underlying degradation process and stochastically changing covariates. The 
present paper falls within such a line of investigation, aiming to include the presence of 
competing risks into the classical shock models scheme. 

In our setting we consider systems characterized by two types of shocks and assume that 
failures are due to a single cause, with shocks arriving according to a pair of independent 
counting processes Nx(t) and N 2 (t). We study in detail the quantities of interest under three 
different failure schemes. In each scheme, failures are assumed to occur at the first instant in 
which one of the following relations holds: 

(a) N 1 (t)+N 2 {t) = M, 

(b) mm{N 1 (t),N 2 (t)} = M, 

(c) max{Ni(t),N 2 (t)} = M, 

where M is a random counting number. We examine various cases arising from suitable 
choices of the distribution of M. 

In Section 2 we describe the general setting of the model by introducing the relevant 
hazard rates and obtain failure densities, survival function, failure probability and the failure 
time moments conditional on cause or type of failure <5. In Section 3 we study this model 
for each of the above schemes under the assumption that Ni(t) and N 2 (t) are homogeneous 
Poisson processes. When failures are due to sum of shocks, failure time T and type of failure 
8 are shown to be independent, and an approximation of the survival function in the presence 
of an arbitrary (large) number h of types of shocks is given. Finally, in Section 4 we suggest 
how to extend this model to the case of two non-exclusive types of shocks. 

2 A model with two kinds of shocks 

Recently, various investigations have been oriented towards multidimensional shock models. 
See, for instance, Ohi and Nishida [27], [2S], and Balu and Sabnis |2J for various results on 
bivariate shock models, including properties of the joint survival probability and preservation 
of reliability structures, and Wong [35] and Pellerey [30] for preservation of stochastic orders 
in multivariate shock models with underlying counting processes. A multidimensional shock 
model driven by counting process is considered in Li and Xu [2T , where each shock may simul- 
taneously destroy a subset of the components of a system consisting of several components. 
We shall now introduce a new kind of shock model characterized by two types of shocks, the 
extension to a higher number of types of shocks being straightforward. This has been inspired 
by the shock model with aftereffects treated in Marshall and Olkin [55] and in Ghurye and 
Marshall [TTj . 

Let T be an absolutely continuous non-negative random variable describing the random 
failure time of a system or of a living organism. We set 8 — i if the failure occurs due to a 
shock of type i, for i = 1,2. Let N(t) denote the total number of shocks occurring in [0,i], 
with t > 0. We have 

N(t) = Ni(t) + N 2 (t), t>0, 

where iVj(i) is the counting process representing the number of shocks of type i occurring in 
[0,t],t = l ) 2. 

The state-space Nq oi{N\(t),N 2 (t)) is partitioned into non-empty subsets Sk, k = 0, 1, 2, . . ., 
that will be called failure sets, Sq including the numbers of non-fatal shocks. In other terms, 
if (Ni(t), N 2 (t)) G So then the shocks that arrived up to time t do not cause a failure. Let 
M be the integer- valued random variable that represents the index of which failure set Sk, 
k = 1,2,..., containes the numbers of fatal shocks. This means that if M = k then the system 
fails at the first instant t > in which (Ni(t), N 2 (t)) £ Sk- The probability distribution and 
the survival probability of M will be respectively denoted by 

(1) p k = P(M = k), fc=l,2,..., 
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and 

(2) P k = P(M>k), k = 0,1,2,.... 

For k = 1, 2, . . . we now define 

§P = {(x 1 ,x 2 ) en 2 -S k : (xi + l,x 2 ) G S k }, 

(3) 

S { k 2) = {{x u x 2 ) e N 2 - 5 fc : (zi,x 2 + 1) e S k }. 

These will be called risky sets, because Sj: \ i — 1,2, containes all states of Ng that lead to 
a failure at time t at the next occurrence of a shock of type i, given that M — k. Note that 
definitions ([3|) assume that failures are due to a single cause. Moreover, note that in general 
S% is different from Sk-i- 

Denoting by frit), t > 0, the probability density function of the failure time T, we have 

(4) f T (t) = fi(t) + h(t), t>0, 
where fi(t) is the sub-density defined by 

fi(t) = ^P{T<t, 5 = i}, t>0, i = l,2. 
Recalling that S = i if the failure occurs due to a shock of type z, we have 

(5) P(S = i)= / fi(t)dt, i = l,2. 

Jo 

In order to express the sub-densities fi(t), i = 1,2, in terms of the joint probability 
distribution of (Ni(t), -^(t)) we now introduce the hazard rates 

r 1 (x 1 ,x 2 ;t) = lim -P{Ni(t + r) = Xi + 1, N 2 (t + r) = x 2 \Ni(t) = x u N 2 (t) = x 2 \, 

(6) 

r 2 (xi,x 2 ;t) = lim -PiN^t + r) = xi,N 2 (t + r) = x 2 + l\Ni(t) = x x ,N 2 {t) = x 2 }, 

with (ari, x 2 ) € Nq and t > 0. Given that x\ shocks of type 1 and x 2 shocks of type 2 occurred 
in [0, t], Ti(x\, x 2 ; t) gives the intensity of the occurrence of a shock of type i immediately 
after t, with i = 1,2. We note that r\(x\, x 2 ; t) + r 2 {x\, x 2 ; t) is the hazard rate of a shock of 
type 1 or 2. 

From the above assumptions it follows that the system fails around time t due to a shock 
of type i if the two-dimensional counting process (Ni(t), N 2 (t)) takes values in at time t 
and a shock of type i occurs immediately after. Hence, conditioning on M and recalling (fTj). 
([3|) and ((6|), for t > and i = 1,2, failure densities can be expressed as 

(7) /i(*) = X>* E P{N 1 (t)=x 1 ,N 2 (t)=x 2 } ri(xi,x 2 ;t). 

A relation similar to (0 holds for the survival function of T, denoted by 

F T {t) = P{T > t}, t>0. 
Indeed, conditioning on (N±(t), N 2 (t)) and recalling ((2]), we obtain 

(8) F T (t) = J2 T k 12 P{Ni{t)=x u N 2 {t)=x 2 }, t>0, 

k=0 (x 1 ,x 2 )€S k 
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where Pq = 1 . Other quantities of interest are the moments of the failure time conditional 
on the cause of failure: 

i r + °° 

(9) E{T S \8 = i) = — / i s /i *)di, * = 1,2. 

p {° = V Jo 

This can be evaluated making use of ((SJ) and 0. 

3 Shocks driven by a Poisson process 

In this Section we assume that the two kinds of shocks affect the system according to inde- 
pendent homogeneous Poisson processes Nx(t) and A^(t). Therefore, for x\,x% — 0, 1,2, . . . 
and t > one has: 

(10) P{7V 1 (t)=x 1 ,iV 2 (t)=.T 2 } = L2J_ 

with Ai > and A 2 > 0. Under assumption (jTDJ) , from ([S]) we now have 

r t (xi,x 2 ;t) — Aj, (an, #2) 6 Nq, £ > 0, £=1,2. 
Hence, from © and © we obtain: 

(U ) m*) = X> £ !^!M!i£±!M!l , i , t >„, 

-f« 2:11 2:21 

and 

(12) e (£0 . 

fc=0 (ii,i 2 )ES|, 

Recalling ([5]) and (fTTj) . it is not difficult to see that 

(13) P(5 = ? :)=7T^ Pfe £ ^if- ^TTf', 1=1,2, 

k=1 (xuX 2 )£Si l) 

where we have set 7Tj = A,/(Ai + A 2 ), i = 1,2. Similarly, from ([9]) we obtain the conditional 
moments 

{xi + x 2 + sy. Xl X2 



j j ^2 

The above expressions show that the structures of the failure sets and of the risky sets are 
essential to specify the nature of the shock model. Hereafter we consider three special cases 
of interest, in which these sets are chosen according to typical models of reliability theory, as 
sketched in Figure Q] for Sk ■ 
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Figure 1: Failure sets Sk when failures are due to (a) sum of shocks, (b) minimum of shocks, 
and (c) maximum of shocks. 

3.1 Failures due to sum of shocks 

Assume that the system fails when the sum of shocks of type 1 and of type 2 reaches a random 
threshold that takes values in {1, 2, . . .}. In this case we have 



S k = {(xi,x 2 ) S Nl : xi + x 2 = k}, k = 0, 1, 2, . 



and 



S 



(1) _ 5(2) _ 



S k — Sk-i, 



fc = l,2, 



Under these assumptions, from (TlTj) and by use of Newton's binomial theorem for i = 1, 2 we 
obtain 



(15) 



/2 (t) =Aie -^+^5:p fe - [(Ai+A2)<]fe " 



k=l 



Similarly, making use of (fT2")) we have 



(16) 



F T (t) 



- „-(Ai+A 2 )t 



fe=0 



(fc-1) 



[(Ai + A 2 )t] A 
kl 



Moreover, from (fT5|) the probability that the failure ultimately occurs due to a shock of type 
i follows: 



(17) 



P(5 = i) 



A; 



Ai + A 2 



1,2. 



Finally, from (fbi|) it is not hard to prove that the conditional moments for s = 1, 2, ... in this 
case are given by 



(18) 



£(T S \S = i) = 



(Ai + A 2 > 



£[(M + s-l) s ], 



1,2, 



where (m) s denotes the descending factorial m(m — l)(m — 2) • • • (m — s + 1). Note that the 
right-hand-side of (fT8|) does not depend on (5. This is not surprising, since the time of failure 
T and the cause of failure S are independent if the system failures are due to sum of shocks, 
as will soon be proved. It is interesting to recall that the relevance of independence between 
T and 5 has been pointed out by various authors (see, for instance, Carriere and Kochar [4]). 



Proposition 3.1 T and S are independent under the model assumptions of Section \3.1\ 
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(ii) 


(fc-1)! 


A i e-"- T / (2 v ^ r ) 




eg; 

n = 


n/2 / \ 


(iii) 


1 


^{l- e -(l + r)} 


1 - e" 


-T 




fe(fc + l) 


7" 






(iv) 


kp 2 (l-p) h ~ 1 


A iP 2 {l + (l-p)r|e-^ 




l+p(l- 


- P) T ) 



Tabic 1: Subdensities and survival function for the model with failures due to sum of shocks, 
with t = (Ai + \2)t. 



Proof. From (|15p and (fTT)) . expressing pk as Pfc-i — Pk there holds: 

P(, ==< ) Mt ) = Ai e-C^)«{g g p J(Ar + A 2 )^ _g g ft ^+^f" 1 } 

^ fe=0 7=fc+l ' k=lj=k+l * >' } 



X -(A 1+ A 2 )t 



Oj'=fc+l fe=l j'= 

„ ^ [(Ax + A 2 )t] fc +^ ^ [ (Al + A 2 )t] fe " 1 



^'2v fc! 2^'^ (fc _i)i 

j=l fc=0 j=2 k=l v ' 



Aie -(A 1+ A 2)t g p /A 1 + A 2 >]i-^ = ^ . = i2 



The proof is thus completed, for T and 5 are independent if and only if /j(t) = P(<$ = ijfrit) 
forallt>0andi = l,2. ■ 

For the model in which failures are due to sum of shocks, the subdensity fi(t) and the 
survival function Fr(t) are shown in Table[T]when the distribution of M is (i) geometric, (ii) 
Poisson over {1, 2, . . .}, (iii) pk — k (k+i) ' ^ = 1» 2, . . ., and (iv) a suitable negative binomial. 
I n (x) denotes the modified Bessel function of the first kind. 

We notice that fi(t) and Pt(£) can be expressed similarly to Eqs. (I15| and (fT6|) if iVi(i) and 
iV2(t) are independent non-homogeneous Poisson processes. However, in this case P(S = i) 
would be equal to the right-hand-side of (|17p only if the ratio of the time-varying intensity 
functions of N\(t) and iV^i) is a constant. Hence, Proposition 13.11 still holds if Ni(t) is a 
Poisson process characterized by an intensity function of the form Ai u(t), with i — 1,2 and 
t > 0. 

The model considered in this paper can be easily extended to the case of an arbitrary num- 
ber h of types of shocks, by introducing a multidimensional counting process (Ni(t), . . . , Nh(t)), 
where Ni(t) describes the number of shocks of i-th type occurred in [0, t], i = 1, 2, . . . , h. In 
this case, under the assumption that failures are due to sum of shocks, the survival function 
of T can be expressed as 

+oo 

(19) F T (t) = J2P kP { Nl ^ + --- + Nh ^ <k }> t ^°- 

k=0 

The form of the right-hand-side of p^|) suggests that a central limit theorem might be imple- 
mented. For instance, assuming that the components of (Ni(t), . . . ,Nh(t)) are independent 
and identically distributed with finite mean and variance, Eq. (|19[) at a fixed time t > for 
large h gives 



where /j,q = E[Ni(to)] and ctq = Var[Ni(to)], with $(■) denoting the standard normal distri- 
bution function. 
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3.2 Failures due to minimum of shocks 



In this Section we consider another failure scheme for the shock model characterized by two 
kinds of shocks. Namely, in this case the system fails when the minimum of the shocks of type 
1 and of type 2 reaches a random threshold that takes values in {1, 2, . . .}. Such a model may 
for instance be appropriate to describe the failure of systems made out of units connected in 
parallel. 

In this model the failure sets and the risky sets are respectively given by 
(20) S k = {(x 1 ,x 2 )enl:mm(x 1 ,x 2 ) = k}, fe = l,2,..., 

and 

= {(* - 1, x 2 ) e : x 2 = k, k + 1, . . .}, 

S { * ] = {(xi,k - 1) G Nl : xi = k, k + 1, . . .}. 
From (fTTj) and (|20[) the following subdensities now follow: 

fi(t) = A ie -(^)* £ Pfc i^L_£ fe (A 2 i), t > 0, 



(21) k=1 

h(t) = \ 2 e-^+^ ^ Pk Wi> E k (X lt ), t > 0, 



k=l 

where for k = 1,2, . . . we have set 



The survival function of T can now be easily obtained from (fT2"|) and (j2"0|) : 

(22) F T (t) = e -(^+A 2 )tgpJM!^ (Alt) + 4 > o. 



Tables [2] and [3] show respectively fi(t) and Fx(t) for the present model, for three choices 
of the distribution of M. Note that the function 



M ntk (x) = 



+°° ^r+n+fc 



— ' r\(r + n)!(r + fc)! ' 



r=0 



appearing in case (ii), is a modified two-index Bessel function (also known as a modified 
Humbert function, see for instance Dattoli et al. [7]). It appears also in the probability dis- 
tributions of certain two-dimensional random motions (see [5] and [5]). The survival function 
in case (ii) of Table [3] has been obtained by recalling the relation (see, for instance, Eq. (5.3) 
of 0) 

+00 + 00 s / 1 \ 1 

J2 J2 an p m M n,M = ex p U ( « + /3 + -3 ) \- 

n=0 m=0 ^ P / J 



3.3 Failures due to maximum of shocks 

Let us now consider a failure scheme in which the system fails when the maximum of the 
shocks of type 1 and of type 2 reaches a random threshold taking values in {1, 2, . . .}. For 
instance, this model is suitable to describe the failure of systems composed by units serially 
interconnected. 
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fi(t) 



(l-p)A; 



7, l+ i(at) 



+°° / 02 \ ™+i 



(Hi) - jJ-^z^a^-^ + E 



A; V / r! (r + 1)! 



Table 2: Subdensities for the model with failures due to minimum of shocks, for cases (i)-(iii) 
of Table 1, with r = (Ai + A 2 )i, a = 2^{l-p)X 1 X 2 and /? = (17A1 A 2 i 2 ) 1/3 . 



F T (t) 



(0 ^ T E 

n=0 



Ai 



(l-p)Aa 



/n(at) + 



A 2 



(l-p)Ai 



(ii) 



cxp < /3 



VA 2 r;y 



Aft^y /^2_ 
AiA 2 tV 1 V V J VAir; 2 

Air; 



+ exp { p 



E^) 



j=0 



* I jtS V 7 (AiA 2 )i +1 V ' VAIA^ 



A2 A2 Ai 



Table 3: Survival function for the model with failures due to minimum of shocks, for the same 
cases of Table 2. 



In this case we assume that the failure sets are given by 
(23) Sk — {(£1,2:2) G No : max(xi,x 2 ) = k}, A: = 1,2,..., 

so that, for k = 1, 2, . . ., the risky sets are 

= {(fc - 1, x 2 ) G N 2 x 2 = 0, 1, . . . k - 1}, 

S<f) = {( Xu fc - 1) g Ng : aji = 0, 1, . . . A - 1}. 

Note that, even if definition (|23|) includes state (k, k) in Sk, the assumption that failures are 
due to a single cause excludes in this model the possibility that states of the type (k, k) are 
failure states. 

From (TTTj) it follows 

+°° (\ fifc-1 

hi t) = x 1 e-(^^Pkjizrw E ^ x ^ 1 ^ °> 

(OA) k=i [ ' >' 

/a(t) = A 2 e-(^)* J^p^-lA—EkiX.t), t > 0, 

where Ek{x), x G R, is the auxiliary function such that J3o(aj) = and 

fe-i j 

(25) £ fc (x) =Y,^\ =e x ~E k {x), fc = l,2,.... 

It is not hard to see that the subdensities fiit) for the two models in which failures are due 
to minimum and maximum of shocks are closely related. This is due to the similar nature 
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of failure sets S k defined in ([20]) and (|23|) . Indeed, making use of identity ([25]) . the sum of 
subdensities (I^T]) and is seen to be 

fe=i ^ >' 

Moreover, from (fT2| and ([23]) the survival function of T is obtained: 

(26) i? T (i) = e-^+^f^pJ^E^t) + ^E k+1 (X 2 t)\, t > 0. 

As for subdensities the survival functions of T for the last two considered models are 

related. Indeed, making again use of ([25]) the sum of functions ([22]) and ([26]) is 

2 +°° / \ +^fe 

E^E^^ *><>. 

i=l fc=o 



4 Concluding remarks 

We have introduced a formulation of shock models driven by bivariate Poisson process that 
includes a competing risks set-up, in which failures are due to one out of two mutually exclusive 
causes. We have obtained specific expressions for the failure densities and the survival function 
under three different failure schemes, in which the variable that causes the failure is (a) the 
sum of shocks, (b) the minimum of shocks, (c) the maximum of shocks. 

Future developments of this model will be oriented to the construction of other structures 
for failure sets Sk, and to the study of a case in which the two fonts of shocks are not mutually 
exclusive, thus including the possibility that failures are due to both types of shocks. The 
underlying counting process will not have independent components; for instance the following 
bivariate Poisson process could be employed (see Marshall and Olkin and references 
therein): 

min(x u x 2 ) ^ k k k +X2 _ k 

P{N 1 (t)=xi,N 2 (t)=x 2 } = e-^^ \ \u 3 \ V u > 

^ (£1 — ky. (x2 — ky.kl 

Furthermore, the structure of the failure sets and of the risky sets will be modified accordingly, 
by taking also into account the necessity that a third hazard rate must be added to those 
defined in Eqs. ([6]). 



Acknowledgments 

The authors thank Petr Lansky, Franco Pellerey, Luigi M. Ricciardi and Hiromi Seno for useful 
comments and fruitful discussions. This work has been performed under partial support by 
MIUR (PRIN 2005), G.N.C.S.-INdAM and Regione Campania. 



References 

[1] (2006) J. Statist. Plann. Inference Special Volume 136, n. 5. 

[2] Balu, M.N. and Sabnis, S.V. (1997) Preservation of certain dependent structures 
under bivariate homogeneous Poisson shock models. Stat. Prob. Lett. 35, 91-100. 

[3] Bedford, T. and Cooke, R. (2001) Probabilistic Risk Analysis: Foundations and 
Methods, Cambridge University Press, Cambridge. 



9 



Carriere, K.C. and Kochar, S.C. (2000) Comparing sub-survival functions in a 
competing risks model. Lifetime Data Anal. 6, 85-97. 

Cesarano, C. and Di Crescenzo, A. (2001) Pseudo-Bcssel functions in the descrip- 
tion of random motions. In: Proceedings of the Workshop on Advanced Special Functions 
and Integration Methods, Mem, 18-23 June 2000, (Eds. G. Dattoli, H.M. Srivastava and 
C. Cesarano), 211-226. Aracnc, Roma. ISBN 88-7999-298-8. 

Crowder, M. (2001) Classical Competing Risks. Chapman & Hall, Boca Raton. 

Dattoli, G., Lorenzutta, S., Maino, G., Torre, A., Voykov, G. and Chiccoli, 
C. (1994) Theory of two-index Bessel functions and applications to physical problems. 
J. Math. Phys. 35, 3636-3649. 

Di Crescenzo, A. (2002) Exact transient analysis of a planar random motion with 
three directions. Stoch. Stoch. Rep. 72, 175-189. 

Di Crescenzo, A. and Longobardi, M. (2006) On the NBU ageing notion within 
the competing risks model. J. Stat. Plann. Infer. 136, 1638-1654. 

El-Neweihi, E., Proschan, F. and Sethuraman, J. (1983) A multivariate new 
better than used class derived from a shock model. Oper. Res. 31, 177-183. 

Esary, J.D., Marshall, A.W. and Proschan, F. (1973) Shock models and wear 
processes. Ann. Prob. 1, 627-649. 

Fagiuoli, E. and Pellerey, F. (1993) New partial orderings and applications. Naval 
Res. Logist. 40, 829-842. 

Fagiuoli, E. and Pellerey, F. (1994) Mean residual life and increasing convex com- 
parison of shock models. Stat. Prob. Lett. 20, 337-345. 

Fagiuoli, E. and Pellerey, F. (1994) Preservation of certain classes of life distribu- 
tions under Poisson shock models. J. Appl. Prob. 31, 458-465. 

Ghosh, M. and Ebrahimi, N. (1982) Shock models leading to increasing failure rate 
and decreasing mean residual life survival. J. Appl. Prob. 19, 158-166. 

Gottlieb, G. (1980) Failure distributions of shock models. J. Appl. Prob. 17, 745-752. 

Ghurye, S.G. and Marshall, A.W. (1984) Shock processes with aftereffects and 
multivariate lack of memory. J. Appl. Prob. 21, 786-801. 

Gut, G. (1990) Cumulative shock models. Adv. Appl. Prob. 22, 504-507. 

Kebir, Y. (1994) Order-preserving shock models. Prob. Engin. Inf. Sci. 8, 125-134. 

Kijima, M. and Nakagawa, T. (1991) A cumulative damage shock model with im- 
perfect preventive maintenance. Naval Res. Logist. 38, 145-156. 

Klefsjo, B. (1981) Survival under the pure birth shock model. J. Appl. Prob. 18, 
554-560. 

Kochar, S.C. (1990) On preservation of some partial orderings under shock models. 
Adv. Appl. Prob. 22, 508-509. 

Lehmann, A. (2006) Dcgradation-threshold-shock models. In Probability, Statistics and 
Modelling in Public Health, 286-298, Springer, New York. 



10 



[24] Li, H. and Xu, S.H. (2001) Stochastic bounds and dependence properties of survival 
times in a multicomponent shock model. J. Multivariate Anal. 76, 63-89. 

[25] Marshall, A.W. and Olkin, I. (1966) A multivariate exponential distribution. J. 
Amer. Statist. Soc. 62, 30-44. 

[26] Marshall, A.W. and Olkin, I. (1967) A generalized bivariate exponential distribu- 
tion. J. Appl. Prob. 4, 291-302. 

[27] Ohi, F. and Nishida, T. (1978) Bivariate shock models and its application to the 
system reliability analysis. Math. Japon. 23, 109-122. 

[28] Ohi, F. and Nishida, T. (1979) Bivariate shock models: NBU and NBUE properties, 
and positively quadrant dependency. J. Oper. Res. Soc. Japan 22, 266-274. 

[29] Pellerey, F. (1994) Shock models with underlying counting process. J. Appl. Prob. 
31, 156-166. 

[30] Pellerey, F. (1999) Stochastic comparisons for multivariate shock models. J. Multiv. 
Anal. 71, 42-55. 

[31] Perez-Ocon, R. and Gamiz-Perez, M.L. (1995) On the HNBUE property in a class 
of correlated cumulative shock models. Adv. Appl. Prob. 27, 1186-1188. 

[32] Perez-Ocon, R. and Gamiz-Perez, M.L. (1996) HNBUE property in a shock model 
with cumulative damage threshold. Comm. Statist. Theory Methods 25, 345-360. 

[33] Shaked, M. and Shanthikumar, J.G. (1991) Shock models with MIFRA time to 
failure distributions. J. Stat. Plann. Infer. 29, 157-169. 

[34] Singh, H. and Jain, K. (1989) Preservation of some partial orderings under Poisson 
model. Adv. Appl. Prob. 21, 713-716. 

[35] Wong, T. (1997) Preservation of multivariate stochastic orders under multivariate Pois- 
son shock models. Adv. Appl. Prob. 34, 1009-1020. 

Antonio Di Crescenzo 

DlPARTIMENTO DI MATEMATICA E INFORMATICA, UNIVERSITA DI SALERNO 

Via Pontc don Mclillo, 84084 Fisciano (SA), Italy 
E-mail: adicrescenzo@unisa.it 



Maria Longobardi 

DlPARTIMENTO DI MATEMATICA E APPLICAZIONI, UNIVERSITA DI NAPOLI FEDERICO II 

Via Cintia, 80126 Napoli, Italy 
E-mail: mar ia . longobardiOunina . it 



11 



